Infrared phonon dynamics of multiferroic BiFeC>3 single crystal 
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We discuss the first infrared reflectivity measurement on a BiFeOa single crystal between 5 K 
and room temperature. The 9 predicted afe-plane E phonon modes are fully and unambiguously 
determined. The frequencies of the 4 A\ c-axis phonons are found. These results settle issues 
between theory and data on ceramics. Our findings show that the softening of the lowest frequency 
E mode is responsible for the temperature dependence of the dielectric constant, indicating that 
the ferroelectric transition in BiFeC>3 is soft-mode driven. 

PACS numbers: 78.30.-j, 77.84.-s, 77.22.-d 
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The interplay between magnetic and dielectric proper- 
ties is an intriguing subject that was already discussed 
by Pierre Curie in 1894 Such, magnetoelectric ef- 
fects were subject of extensive studies during the 1960s 
@; 0, Sj- More recently, the ability of controlling the 
dielectric properties through a magnetic field and vice- 
versa renewed the interest in these materials @, 0, S|- 
Coupling between electrical and magnetic properties are 
particularly interesting in multiferroic materials, which 
present simultaneously two or more ferroic or antiferroic 
order parameters. The ferroelectric Q and antiferro- 
magnetic properties of BiFeOa based systems, showing 

T Curie = 1090 K @ and T Neel = 673 K 13' have been 
known for a long time but the possibility of making de- 
vices operating at room temperature [9J] gave a new burst 
to bismuth ferrite. The coexistence of ferroelectricity and 
antiferromagnetism at room temperature, makes BiFeOa 
one of the most interesting multiferroic materials. 

Optical spectroscopy is a very powerful tool to under- 
stand the driving mechanism of the ferroelectric transi- 
tion and, eventually, its coupling to ma gne tic ordering 
ToL 11 1 . Nevertheless, only a few Raman 12, 7J|, Ta\ and 
infrared (IR) [TJI, [Tj| studies on BiFeOa are known to 
date. Surprisingly none of the IR measurements were car- 
ried out on single crystals. There are two important open 
questions concerning the electrodynamics of BiFeOa: (i) 
the phonon modes symmetry; and (ii) the role (or ab- 
sence of thereof) of phonon softening in the ferroelectric 
transition. 

In its ferroelectric phase, BiFeOa belongs to the R3c 
(C® v ) space group, which derives from the ideal (para- 
electric) Pmirn (OjJ cubic perovskite group by a small 
distortion along the (111) cubic directions. Its magnetic 
structure is a spiral-cycloidal incommensurate phase. 
Group theory analysis for the R3c ferroelectric group pre- 
dicts 4Ai ® 5A 2 © 9E optical phonon modes [Tflj. The A\ 
modes are infrared active along the c-axis, the E modes 



are infrared active on the afe-plane and the A-i modes are 
silent. Early Raman data [12j found at least 11 modes 
at room temperature. Kamba et al. [l6| proposed to 
fit the infrared data on ceramics using 13 Lorentz os- 
cillators. Although these results found the proper num- 
ber of modes, they cannot clarify the symmetry of the 
phonons. Some extra insight on the phonon assignment 
was obtained on more recent Raman data 0, Hi] but 
some inconsistencies in phonon assignment remain. In 
addition, Haumont et al. 12] argued that because the fre- 
quency of the lowest mode did not vanish smoothly at the 
Curie temperature, the ferroelectric transition in BiFeOa 
would not be soft-mode driven. This picture was later 
revised to propose that the incomplete phonon softening 
could rather be a sign of a first order transition 1|| 2(| . 
However, both scenarios remain speculative. 

In this paper, we show the polarized infrared response 
from 5 K to room temperature of a BiFeOa single crys- 
tal. Our crystal geometry allows us to unambiguously 
determine the afe-plane E phonons and to infer the fre- 
quencies of the c-axis A\ phonons. Moreover, we show 
that the temperature variation of the dc dielectric func- 
tion is completely governed by the softening of the lowest 
frequency mode giving a strong evidence for a soft-mode 
driven transition. 

High quality and purity BiFeOa crystals were obtained 
by flux growth 17|,ll8(. We measured a thin BiFeOa plate 
measuring 0.8 x 0.8 x 0.01 mm 3 . Near normal incidence 
infrared reflectivity was recorded from 25 to 4000 cm -1 
in a Bruker IFS66v interferometer. Below 650 cm -1 we 
utilized a Ge coated Mylar beamsplitter and a 4 K Si 
bolometer. Above 500 cm -1 a Ge:KBr beamsplitter and 
a MCT detector were used. In the overlapping region, the 
spectra agreed within 1 %. The sample was mounted on 
the cold finger of an ARS Helitran open cycle cryostat 
equipped with polyethylene and KRS-5 windows. The 
spectra were collected on the as-grown 0.8 x 0.8 mm 2 flat 
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FIG. 1: (Color online) In-plane (electric field of light E _L z) 
infrared reflectivity of a BiFe03 single crystal at 5 K and 300 
K. The dashed lines are the data and the solid lines are fits 
using a dielectric function in the form of Eq. [1] The lower 
arrows indicate the TO frequencies used in our fits. The upper 
arrows are the frequencies calculated in Ref. [2(1 

face of the sample. To obtain the infrared reflectivity on 
such a small sample, an overfill technique with a in-situ 
gold evaporation was performed |19( . The measured sur- 
face, corresponding to the (012)h cx plane, is reminiscent 
of the high temperature cubic phase [(010) cu b]- X-ray 
analysis on similar crystals, showed that the axis per- 
pendicular to the sample flat surface makes an angle of 
~ 55° with the c-axis [17| ■ Using polarized light we could 
isolate the contribution from the afe-plane. Our crystal 
geometry does not allow to isolate the c-axis response 
which is heavily contaminated by the a&-plane phonons. 

Figure [T] shows the infrared reflectivity at 5 K and at 
room temperature for the electric field polarization lying 
in the a&-plane. The dashed lines are the data and the 
solid lines are fits using a dielectric function of the form 

£ (U) = Eoo + > -2 2 . , (1) 

where ftroj 7 lj an d are , respectively, the transverse 
optical (TO) frequency, damping and oscillator strength 
for the j-th mode. £00 is the permittivity due to high 
frequency electronic excitations. The near normal inci- 
dence reflectivity is obtained from the Fresnel formula 
R = \(V^ ~ + -01 ■ The fitting parameters for 

these two temperatures are shown in Table HI 

The 5 K spectrum in Fig. Q] clearly shows the 9 E 
modes predicted by group theory. At room tempera- 
ture, two of these phonons [E(3) and E(5)] are strongly 
damped but keeping them improves the fit quality. Her- 
met et al. [20| calculated the phonon frequencies ex- 
pected for BiFeOs using density functional theory. The 



TABLE I: Fitting parameters used in Eq. \T\to describe the in- 
plane (E _L z) infrared reflectivity of a BiFeOs single crystal. 
The high frequency dielectric function (Sao) was kept constant 
at 9.04 for all temperatures. Frequencies and dampings are 
in cm -1 and the oscillator strengths are dimensionless. 
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lower arrows (black) at the bottom of Fig. [T] show the 
TO frequencies obtained from our fit at 5 K. The up- 
per arrows (green) indicate the values calculated in Ref. 
I20L For most modes, the agreement is excellent (better 
than 5 %). Modes E(l) and E(2), however, are found 
at frequencies 20-30% smaller than the calculated val- 
ues. Conversingly, mode E(8) is found at an energy 8 % 
higher than the theoretical value. Compared to the data 
reported by Kamba et al. in a ceramic material fl6| , our 
obtained frequencies are slightly upshifted, but the most 
important result concerns the modes assignment, once 
we only agree in the case of five modes, E(l), E(4), E(6), 
E(7) and E(9). The use of polarized light allowed us to 
correctly assign the other modes. 

Figure [2] shows the temperature evolution of the TO 
frequencies for the 9 E phonon modes. The E(l) mode 
softens strongly from 5 K to room temperature. High 
temperature infrared [l6| and Raman [l2 . HH measure- 
ments indicate that this mode continues to soften but 
instead of smoothly decreasing to zero frequency it drops 
rather abruptly at the transition. Modes E(2), E(6) and 
E(8) show the same behavior as mode E(l). Interest- 
ingly, mode E(7) hardens with increasing temperature. 
The frequencies of the remaining modes have little or no 
temperature dependence. According to the calculated 
phonon frequencies j^, E(l) and E(2) are related to vi- 
brations involving the Bi atoms. E(6) and E(8) should be 
dominated by oxygen motion. Because of their softening, 
it is reasonable to think that these modes also involve mo- 
tion of Bi atoms. Concerning E(7) — also dominated by 
oxygen motion — its hardening with increasing temper- 
ature shows that this mode has a particular dependence 
on the structural evolution, i.e., it is more sensitive to 
the decreasing rhombohedral elastic distortion than to 
the unit cell volume expansion, upon heating. 

More insight on the E(l) mode can be obtained from 
an analysis of its oscillator strength. One can calculate 
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FIG. 2: Temperature dependence of the in-plane (E _L z) 
transverse optical frequencies for a BiFeOa single crystal. Due 
to its very small oscillator strength, the parameters for mode 
E(3) have been kept constant above 200 K. Unless indicated 
otherwise, the values uncertainties are within the symbols 
size. 



the dielectric constant from the optical contributions as 



6(0) 



3 = 1 



(2) 



This calculation, for our results, is shown as open trian- 
gles in Fig. [3l Our values are twice the values obtained by 
Kamba et al. [16[ . This is a consequence of the higher re- 
flectivity and absence of microstructure (porosity, grain 
boundaries) in our sample. The values in Fig. [3] are 
a bit smaller (10 %) than those obtained by dielectric 
measurements at 1 GHz on ceramics solid solutions of 
BiFe0 3 with Pb(Ti,Zr)0 3 0. We also plotted in this 
figure the temperature dependence of the E(l) mode os- 
cillator strength (solid circles). It becomes clear that the 
temperature increase in e(0) is explained by the softening 
of the E(l) mode alone. Indeed the /-sum rule states that 
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FIG. 3: (Color online) The open triangles are the temperature 
dependence of the dielectric constant calculated using Eq. [2] 
(left axis). The solid circles depict the thermal evolution of 
the oscillator strength for the E(l) mode (right axis). As both 
axis have the same range for e, one can see that the tempera- 
ture dependence of the dielectric constant comes mostly from 
the softening (hence oscillator strength increase) of the E(l) 
phonon. 



the sum of plasma frequencies squared Ae,-^^ -) 
should be a constant independent of temperature [2 if . 
An even more stringent condition is that, in the ab- 
sence of phonon coupling, each individual plasma fre- 
quency should remain constant. Indeed, within ±5%, the 
plasma frequency for the mode E(l) respects this condi- 
tion. Therefore, the softening of the E(l) phonon induces 
an increase in its oscillator strength which can fully ac- 
count for the increase in the dielectric constant. This is 
a strong indication that this phonon participates actively 
as the soft mode driving the ferroelectric transition. As 
this mode dies sharply at the Curie temperature instead 
of vanishing smoothly [12, it was first considered 
that the transition was not soft-mode driven [12]. This 
picture was later revised to propose that the transition 
was displacive and did have a soft mode which incom- 
plete softening could be a sign of a first order transtion 
[16L [20I] . Figure 3, gives a strong support for this latter 
proposition. 

An important point arises when we look at the theo- 
retically obtained oscillator strengths and our experimen- 
tally determined values. Hermet et al. [2(| calculated the 
phonons TO and longitudinal optical (LO) frequencies 
and found SItoi — 102 cm" 1 and Qloi = 104 cm -1 for 
mode E(l) and flroi = 152 cm -1 and SIloi = 175 cm" 1 
for mode E(2). We can use the TO-LO splitting to cal- 
culate each mode oscillator strength [22]. The TO-LO 
splittings calculated by Hermet et at, yields theoretical 
oscillator strengths of Aei = 2.1 and AE2 = 12.3. This 
is strikingly different from the values we determined ex- 
perimentally (see Table J]). We find that the oscillator 
strengths for the E(l) and E(2) modes are of the same 
order of magnitude and, more important, the oscillator 
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FIG. 4: (Color online) The thinner line is the in-plane (E _L z) 
reflectivity. The thicker line is the combined response. The 
upward pointing arrows are the TO frequencies calculated in 
Ref. [U The down pointing arrows are our proposition for 
the E || c phonons. 

strenght of mode E{1) is large enough (by a factor 10 
with respect to the calculation) to drive the changes in 
the dielectric constant. Our results, opposite to Hermet 
et al. and in agreement with Kamba et ai, show that the 
E(l) mode is the soft mode in BiFeOa. 

Although we cannot measure directly the c-axis re- 
sponse we can choose a polarization direction in our mea- 
surement crystal face that minimizes the a&-plane con- 
tribution. We can thus infer the phonon locations along 
the c-axis. 

In Fig. 5]we compare the afe-plane response (thick gray 
line) to the mixed spectrum (thin red line) obtained at 
5 K, temperature at which the phonons are sharp and 
well defined. There are three phonon modes that are 
very clearly present in the mixed response and absent 
from the a6-plane spectrum. They are indicated by the 
down pointing arrows. The forth mode is very likely ap- 
pearing as a broader shoulder just above the E{2) mode 
(~ 180 cm -1 ) indicated by the arrow with a question 
mark. We assign these features to the 4 A\ modes. The 
upward pointing arrows at the bottom of the figure in- 
dicate the frequencies for the A\ modes calculated from 
first principles [2(| . The first two A\ phonons found in 
our spectra (180 cm -1 and 215 cm -1 ) agree relatively 
well with the predicted values. However, the phonon ex- 
pected around 320 cm -1 is found closer to 475 cm -1 . The 
phonon expected at 520 cm -1 appears above 600 cm -1 . 
The disagreement for the latter phonon can be under- 
stood in terms of the strong E(9) phonon contamination 
in our mixed response spectra. Because of this large E(9) 
mode, the Ai (4) mode may appear as a dip at its LO fre- 
quency, rather than a peak close the TO energy. The cal- 
culated LO frequency for the ^4i(4) mode is ~ 535 cm^ 1 , 



closer to our finding but still a poorer agreement than the 
first two modes. We also remark that the assignment of 
the Ai(3) and Ai(4) modes differ substantially from pre- 
vious propositions ITJ. Pal , although A\(3) agrees with 
recent Raman results [14j |. 

We measured the low temperature polarized infrared 
response of a BiFeOa single crystal. We determined 
clearly all the afr-plane phonon parameters. Our find- 
ings agree very well with calculated phonon frequencies 
and solve the phonon symmetry assignment issues. The 
c-axis phonons could not be completely determined from 
our crystal geometry but their frequencies are inferred 
from a mixed polarization measurement. We show that 
the softening of the E{\) phonon explains completely the 
temperature dependence of the dielectric constant. Our 
results are a clear evidence that the E(l) phonon is un- 
doubtedly the soft mode driving the ferroelectric transi- 
tion. 
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